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ABSTRACT 

We investigate torsional Alfven oscillations of relativistic stars with a global dipole magnetic 
field, via two-dimensional numerical simulations. We find that a) there exist two families of 
quasi-periodic oscillations (QPOs) with harmonics at integer multiples of the fundamental 
frequency, b) the lower-frequency QPO is related to the region of closed field lines, near the 
equator, while the higher-frequency QPO is generated near the magnetic axis, c) the QPOs are 
long-lived, d) for the chosen form of dipolar magnetic field, the frequency ratio of the lower to 
upper fundamental QPOs is ~ 0.6, independent of the equilibrium model or of the strength of 
the magnetic field, and e) within a representative sample of equations of state and of various 
magnetar masses, the Alfven QPO frequencies are given by accurate empirical relations that 
depend only on the compactness of the star and on the magnetic field strength. The lower and 
upper QPOs can be interpreted as corresponding to the edges or turning points of an Alfven 
continuum, according to the model proposed by Levin (2007). Several of the low-frequency 
QPOs observed in the X-ray tail of SGR 1806-20 can readily be identified with the Alfven 
QPOs we compute. In particular, one could identify the 18Hz and 30Hz observed frequencies 
with the fundamental lower and upper QPOs, correspondingly, while the observed frequencies 
of 92Hz and 150Hz are then integer multiples of the fundamental upper QPO frequency (three 
times and five times, correspondingly). With this identification, we obtain an upper limit on the 
strength of magnetic field of SGR 1806-20 (if is dominated by a dipolar component) between 
~ 3 and 7 x IG^^G. Furthermore, we show that an identification of the observed frequency 
of 26Hz with the frequency of the fundamental torsional £ — 2 oscillation of the magnetar's 
crust is compatible with a magnetar mass of about 1.4 to 1.6AfQ and an EOS that is very stiff 
(if the magnetic field strength is near its upper limit) or moderately stiff (for lower values of 
the magnetic field). 

Key words: relativity - MHD - stars: neutron - stars: oscillations - stars: magnetic fields - 
gamma rays: theory 



1 INTRODUCTION 

The phenomenon of Soft Gamma Repeaters (SGRs) may allow us 
in the near future to determine fundamental properties of strongly 
magnetized, compact stars. Already, there exist at least two sources 
in which quasi-periodic oscillations (QPOs) have be en observed 
in the ir X-ray tail, following the initial discovery by llsrael et aEI 
( |2005|) . see|Watts & Strohmaver (2006) for a recent review. The 
frequency of many of these oscillations is similar to what one 
would expect for torsional modes of the solid crust of a com- 
pact star. This observation is in support of the proposal that SGRs 
are magnetars (compact objects with very strong magnetic fields) 
dOuncan & Thom pson 1992). During an SGR event, torsi onal os- 
cillations in the solid crust of the star could be excited touncani 
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1 19981) . leading to the observed frequencies in the X-ray tail. How- 
ever, not all of the observed frequencies fit the above picture. For 
example, the three lowest observed frequencies for SGR 1806-20 
are 18, 26, 30Hz. Only one of these could be the fundamental, 
^ = 2, m = torsional frequ ency o f the crust, as the first overtone 
has a much higher frequency. iLevinI ( [2006.) stressed the importance 
of crust-core coupling by a global ma gnetic field and of t he ex- 
istence of an Alfven continuum, while IGlampedakis et al.l ( |2006|) 
considered model with simplified geometry, in which Alfven oscil- 
lations form a discrete spectrum of normal modes, th at could be as - 
sociated with the observed low-frequency QPOS. In lLevinl ( l2007h . 
the existence of a continuum was stressed further and it was shown 
that the edges or turning points of the continuum can yield long- 
lived QPOs. In addition, numerical simulations showed that drifting 
QPOs within the continuum become amplified near the frequencies 
of the crustal normal modes. Within this model. Levin suggested 
a likely identification of the 18Hz QPO in SGR 1806-20 with the 
lowest frequency of the MHD continuum or its first overtone. The 
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above results were obtained in toy models with simplified geometry 
and Newtonian gravity. 

In this Letter, we perform two-dimensional numerical simu- 
lations of linearized Alfven oscillations in magnetars. Our model 
improves on the previously considered toy models in various ways: 
relativistic gravity is assumed, various realistic equations of state 
(EOS) are considered and a consistent dipolar magnetic field is con- 
structed. We do not consider the presence of a solid crust, but only 
examine the response of the ideal magnetofluid to a chosen initial 
perturbation. 

Spherical stars have generally two type of oscillations, 
spheroidal with polar parity and toroidal with axial parity. The ob- 
served QPOs in SGR X-ray tails may originate from toroidal oscil- 
lations, since these could be excited more easily than poloidal os- 
cillations, because they do not involve density variations. In New- 
tonian theory, there have been several investigations of torsiona l 
oscillations in the crust region of neutron stars (see e.g.. lLeel ( l2007l) 
for reference). On the other hand, only few studies have taken 
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general relativity into account (Messios et al. 2001; Sotani et al. 
[2007 a,b; Samuelsson & Andersson 2007; Vavoulidis et al. 2007). 

SGRs produce giant flares with peak luminosities of 10''* - 
10*^ erg/s, which display a decaying tail for several hundred sec- 
onds. Up to now, three giant flares have been detected, SGR 0526- 
66 in 1979, SGR 1900-^14 in 1998, and SGR 1806-20 in 2004. The 
timing analysis of the latter two events revealed several QPOs in 
the decaying tail, whose frequencies are approximately 18, 26, 30, 
92, 150, 625, and 1840 Hz for SG R 1806-20, and 28, 53, 84, and 
155 Hz for SGR 1900-F14, see Watts & Strohm aver (2006). 

In ISotanietalJ ( 120073) (hereafter Paper I), it was suggested 
that some of the observational data of SGRs could agree with the 
crustal torsional oscillations, if, e.g., frequencies lower than 155 
Hz are identified with the fundamental oscillations of different har- 
monic index I, while higher frequencies are identified with over- 
tones. However, in Paper I and above, it will be quite challenging 
to identify all observed QPO frequencies with only crustal torsional 
oscillations. For example, it is difficult to explain all of the frequen- 
cies of 18, 26 and 30 Hz for SGR 1806-20 with crustal models, 
because the actual spacing of torsional oscillations of the crust is 
larger than the difference between these two frequencies. Similarly, 
the spacing between the 625Hz and a possible 720Hz QPO in SGR 
1806-20 may be too small to be explained by consecutive overtones 
of crustal torsional oscillations. 

One can notice, however, that the frequencies of 30, 92 and 
150 Hz in SGR 1806-20 are in near integer ratios. As we will show 
below, the numerical results presented in this Letter are compatible 
with this observation, as we find two families of QPOs (correspond- 
ing to the edges or turning points of a continuum) with harmonics 
at near integer multiples. Furthermore, our results are compatible 
with the ratio of 0.6 between the 18 and 30Hz frequencies, if these 
are identified, as we suggest, with the edges (or turning points) of 
the Alfven continuum. With this identification, we can set an up- 
per limit to the dipole magnetic field of ~ 3 to 7 x lO'^G. If the 
drifting QPOs of the continuum are amplified at the fundamental 
frequency of the crust, and the latter is assumed to be the observed 
26Hz for SGR 1806-20, then our results are compatible with a mag- 
netar mass of about 1.4 to 1.6Mq and an EOS that is very stiff (if 
the magnetic field strength is near its upper limit) or moderately 
stiff (for lower values of the magnetic field). 

Unless otherwise noted, we adopt units of c = G = 1, where 
c and G denote the speed of light and the gravitational constant, 
respectively, while the metric signature is ( — ,+ ,+ ,+). 



2 NUMERICAL SETUP 

The general-relativistic equilibrium stellar model is assumed to be 
spherically symmetric and static, i.e. a solution of the well-known 
TOV equations for and perfect fluid and metric described the line 
element 
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We neglect the influence of the magnetic field on the structure of 
the star, since the magnetic field energy, £m, is orders of magni- 
tudes smaller than the gravitational binding energy, Eg, for mag- 
netic field strengths considered realistic for magnetars, £jv\ /Eg ~ 
10~''(i3/(10'^®G))^. For simplicity, we assume that the magnetic 
field is a p ure dipole (toroidal m agnetic fields will be treated else- 
where, see ISotani et al.l ( l2007ct) . Details on the numerical method 
for constructing the magnetic field, as well as representative figure 
of the magnetic field lines, can be found in Paper I. 

MHD oscillations of the above equilibrium model are de- 
scribed by the linearized equations of motion and the magnetic 
induction equations, presented in detail in Papers I and 11 (we ne- 
glect perturbations in the spacetime metric, as these couple weakly 
to toroidal modes in a spherically symmetric background). The 
perturbative equations presented in Papers I and II can readily 
be converted from an eigenvalue problem to the form of a two- 
dimensional time-evolution problem, by defining a displacement 
y{t, r, 6) due to the toroidal motion (the coefficient of shear vis- 
cosity is set to = 0, as we neglect the presence of a solid crust in 
the present work). The contravariant coordinate component of the 
perturbed four- velocity, Su"^ , is then related to time derivative of y 
through 

Sul" ^e'''dty{t,r,e). (2) 
The two-dimensional evolution equation for y(t, r, 9) is 
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where Att, A20, An, A02, Aio, and ^01 are functions of r and 9, 
given by 
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and where ai (r) and ji (r) are the radial components of the elec- 
tromagnetic four-potential and the four-current, respectively. In the 
above equations, a prime denotes a partial derivative with respect 
to the radial coordinate. 

In our numerical scheme, we employ the 2nd-order, iterative 
Crank-Nicholson scheme. A numerical instability, which sets in af- 
ter many oscillations, was treated by adding a 4th-order Kreiss- 
Oliger dissipation term dKreiss & 01igej[l973h . shown as eo'Diy 
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Figure 1. FFT of the MHD oscillations at S ~ (magnetic axis). Three 
lines in each figure correspond to different radial positions r ~ 0, R/2, 
and R). A fundamental (n = 0) QPO and several overtones (nearly integer 
multiples) are clearly present. 
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Figure 2. Same as Fig. [T] but aX. 9 7r/2 (magnetic equator). A second 
family of QPOs is present. Arrows indicate several continuous parts, of 
which only the first is distinct from the others, which partially overlap. 



in Equation ^ above. We experimented with various values of the 
dissipation coefficient e d and found the evolution to be stable for 
values as small as a few times 10~^. We verified that in this limit, 
the solution becomes independent of the strength of the numeri- 
cal dissipation, as the 2?4 Kreiss-Oliger dissipation operator intro- 
duces an error of higher order than the 2nd-order iterative Crank- 
Nicholson scheme. The numerical grid we use is equidistant, cov- 
ering only the interior of the star, with (typically) 50 radial zones 
and 40 angular zones (we also compared our results to simulations 
with 100x80 points). The boundary conditions are: (a) 3^ = at 
r = (regularity), (b) y^r = at r = i? (vanishing traction) 
(c) 3^,e = at 6* = (axisymmetry) and (d) 3^ = at 6) = 7r/2 
(equatorial plane symmetry of the 1 = 2 initial data). We obtain the 
same results if you use a grid that extends to ^ = tt. We have also 
evolved initial data with ^ = 3, for which the appropriate boundary 
condition at ^ = n /2 is y^e = 0. 



3 ALFVEN QPOS 

As a fist equilibrium model, we examined a polytropic stellar model 
with r = 2, whose mass and radius are M = IAMq and 
R = 14.16 km, respectively, with a dipole magnetic field strength 
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Figure 3. Distribution of effective amplitude of several Alfven QPOs (see 
text for details). The grayscale map varies from white (zero amplitude) to 
black (maximum amplitude). 



of B = Bp = 4 X 10^^ G, where Bp is a typical strength. As ini- 
tial data, we use the numerical eigenfunction of the ^ = 2, m = 
fundamental mode, for the truncated system presented in Paper II. 
The grid size is 100 x 80 and ed = 10~^, while the total sim- 
ulation time is 2s. By computing the FFT of dty (which is pro- 
portional to (5ii'*) at various points inside the star, we obtain the 
following results: Examining the FFT at three different radial loca- 
tions (r ~ 0, R/2, J?) at 6 ~ (magnetic axis), Fig.[T] we observe 
a number of narrow frequency peaks. The strength of the peaks is 
several orders of magnitude larger than the FFT continuum. The 
overtones (n > 0) are nearly integer multiples of the fundamen- 
tal frequency (n = 0). The corresponding FFTs at the same radial 
points, but at 61 ~ 7r/2 (magnetic equator) are shown in Fig. |2] 
where, in addition to the family of frequency peaks observed in 
Fig. [T] one can also observe a second family of frequency peaks, 
mainly at r ~ i?, for which the fundamental frequency is at a ra- 
tio of 0.6 with respect to the fundamental frequency in Fig.[T] The 
fundamental frequency of this family, as well as the first overtone 
(again an integer multiple) are clearly visible, while the other over- 
tones seem to be buried inside a continuum of other frequencies. 

Using Levin's toy model for the Alfven continuum in magne- 
tars, the above numerical results can be interpreted as follows: The 
fundamental frequency peaks of the two families are QPOs gener- 
ated by the edges or turning points of an Alfven continuum (hence- 
forth we call these two frequencies as the fundamental lower and 
upper QPOs and denote them as Lq and Uo). We denote the over- 
tones (which are nearly integer multiples) as L„ and t/„, where 
n > 0. For the chosen form of the magnetic field, the extent of 
the first continuum, between the fundamental Lo and Uo QPOs is 
at lower frequencies than the continuous frequency intervals corre- 
sponding to the overtones, which all overlap partially. 

Because of the Alfven continuum, the phase of the oscilla- 
tions at the QPO frequencies is not constant, but varies throughout 
the star. Since the magnetic field is axisymmetric, the axis is both 
a turning point and edge of the continuum. Near the magnetic axis, 
the background magnetic field varies slowly, which allows for the 
upper QPOs to remain nearly coherent for a large number of os- 
cillations. Since <5m'* is only a coordinate component, we also ob- 
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Figure 4. Time evolution of dtV at the location inside the star where the 
effective amplitude of the fundamental upper QPO. Uq attains its maximum 
value. The amplitude of the oscillations barely diminish with time 



tained FFTs of r sin 6dty (which is proportional to the physical ve- 
locity in a unit basis) at all grid points. In two-dimensional simula- 
tions of fluid modes in non-magnetized stars, it has been shown that 
the amplitude of the FFT of physical variables at every point in the 
numerical grid, at a chosen normal mode frequen cy, is correlated to 
the sh ape of the eigenfunction of a given mode jStergioulas et alj 
l2004h . Similarly, we use the magnitude of the FFT of r sin 6dty 
at every point in the grid to define an effective amplitude for each 
QPO (which, of course, will be time- varying, due to the absorption 
of individual modes in the continuum). In Fig. (3] we show the ef- 
fective amplitude of the fundamental upper QPO, Uq and its first 
overtone Ui as well as for the corresponding lower QPOs Lq and 
Li, obtained after a simulation with a 50 x 40 grid, e_D = 2 x 10~^ 
and a duration of 2. Is. 

The effective amplitude for Uo has a maximum near the mag- 
netic axis, while there are also nodal lines along certain magnetic 
field lines. For the overtone, Ui, there exists an additional nodal 
line, starting perpendicular near the magnetic axis, which divides 
the region of maximum amplitude near the magnetic axis roughly 
in half. Similarly (not shown here), each successive overtone cor- 
responds to an additional "horizontal" nodal line, dividing the re- 
gion of maximum amplitude near the magnetic axis into roughly 
equidistant parts. This agrees well with the fact that the frequencies 
of the overtones are nearly integer multiples of the fundamental 
frequency. 

The effective amplitude of the fundamental lower QPO, Lq, 
is practically limited to a region that is only somewhat larger than 
the region of closed magnetic field lines inside the star. For the first 
overtone Li, a nodal line divides the region of maximum amplitude 
into two parts. 

In Fig.|4]we show the evolution of dty at the location inside 
the star where the effective amplitude of the fundamental upper 
QPO, Uo attains its maximum value. It is evident that the QPO is 
long-lived, since the amplitude of the oscillations barely diminishes 
with time. 



4 REALISTIC EOS AND EMPIRICAL RELATIONS 

We have obtained the lower and upper Alfven QPO frequencies 
for a representative sample of magnetar models with realistic (tab- 
ulated) EOSs and various masses. These equilibrium models and 
their detailed properties have already been presented in Paper I. In 
Table [T] we summarize our numerical results for B — B^, for a 
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Figure 5. Quadratic fits in terms of the compactness of the star, M/ R, of the 
lower and upper fundamental Alfven QPO frequencies, obtained for a rep- 
resentative sample of equilibrium models with various EOSs and masses. 
The magnetic field was set to _B = . 



somewhat smaller sample of models, than the one considered in 
Paper I. As in the case of the polytropic model of Sec.|3] the over- 
tones are nearly integer multiples of the fundamental frequency 
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(10) 
(11) 



with an accuracy on the order of 1% or better. In addition, the ratio 
of the lower to upper QPO frequencies (also shown in Table \Q 
roughly agrees with the value of 0.6 for the polytropic model. A 
quadratic fit in terms of M/R gives: 



fUo 



0.62 



1.08 



+ 4.52 



(12) 



with an accuracy of better than 2%. Similarly, a quadratic fit for the 
ratio of the lower to upper first overtone frequencies yields 
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with a similar accuracy. 

We find that the frequencies of the Alfven QPOs scale linearly 
with the strength of the magnetic field (with an accuracy on the 
order of 1% or better) at least for magnetic field strengths on the 
order of B^. 

Within the representative sample of the equilibrium models 
we display in Table[T] the frequencies of all lower and upper Alfven 
QPOs (with n = 0, 1, 2, ...) are given by the following two empir- 



ical relations 
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with an accuracy of less than about 4%. The quadratic fits in terms 
of M/R, for B = Bfj, is shown in Fig. [5] We emphasize that in all 
of the above relations, the ratio M/R is dimensionless, in gravita- 
tional units (c — G — 1). To restore units, it has to be replaced by 

GM/Rc^. 



Alfven QPOs in Magnetars 5 



Table 1. Frequencies of lower and upper Alfven QPOs and their ratios, for a representative sample of equilibrium models, constructed with various EOSs and 
masses and for a magnetic field strength of B = (see text for details). 



Model M/R (Hz) fu,, (Hz) ratio Jl, (Hz) fu, (Hz) ratio fu^ (Hz) 



A+DH14 


0.218 


15.4 


25.0 


0.616 


30.7 


49.4 


0.621 


74.4 


A+DH16 


0.264 


11.7 


18.3 


0.639 


23.5 


35.7 


0.658 


54.0 


WFF3+DHi4 


0.191 


17.9 


29.8 


0.601 


36.2 


59.2 


0.611 


89.8 


WFF3+DH18 


0.265 


11.7 


18.0 


0.650 


23.5 


35.5 


0.662 


53.3 


APR+DH14 


0.171 


20.4 


34.1 


0.598 


41.3 


68.6 


0.602 


104.6 


APR+DH20 


0.248 


12.8 


20.6 


0.621 


26.0 


40.3 


0.645 


61.0 


L+DH14 


0.141 


23.7 


40.8 


0.581 


47.5 


81.6 


0.582 


123.8 


L+DH20 


0.199 


16.4 


27.8 


0.590 


33.1 


54.7 


0.605 


82.6 



5 CONSTRAINTS ON MAGNETIC FIELD STRENGTH 
AND EOS 

Several of the low-frequency QPOs observed in the X-ray tail of 
SGR 1806-20 can readily be identified with the Alfven QPOs we 
compute. In particular, one could identify the 18Hz and 30Hz ob- 
served frequencies with the fundamental lower and upper QPOs, 
correspondingly, while the observed frequencies of 92Hz and 
150Hz would then be integer multiples of the fundamental upper 
QPO frequency (three times and five times, correspondingly). With 
this identification, our empirical relations Eqs. l ll4t and M5\ con- 
strain the magnetic field strength of SGR 1806-20 (if is dominated 
by a dipolar component) to be between 3 x lO^^G and 7 x lO^^G. 
Furthermore, an identification of the observed frequency of 26Hz 
with the frequency of the fundamental torsional 1 = 2 oscillation of 
the magnetar's crust (Eq. (79) of Paper I) implies a very stiff equa- 
tion of state and a mass of about 1 .4 to 1 .6Mq . For example, for the 
IAMq model constructed with EOS L-l-DH, one obtains the fol- 
lowing frequencies: 2to = 25.8Hz, /ig = 17.5Hz, fu^ = 30.0Hz, 
fui = 90.1Hz and fu^ = 150.2Hz, for B = 2.94 x 10^^ G. 

Alternatively, one could also identify the 18Hz and 30Hz ob- 
served frequencies with overtones (which are also at a near 0.6 ra- 
tio). In this case, the strength of the magnetic field derived above 
is only an upper limit and the actual magnetic field may be weaker 
Then, if one assumes that the observed frequency of 26Hz is due to 
the fundamental 1 = 2 crust mode for a weak magnetic field, our 
numerical data agree best with a 1 AMq model constructed with an 
EOS of moderate stiffness. For example, for the 1.4A/q model con- 
structed with the APR-l-DH EOS one obtains 2*0 = 25.9Hz, = 
17.7Hz, fu^ = 30.0Hz, /u^ = 90.1Hz and fu^ = 150.1Hz, for 
B = 1.77 X 10^^ G. 



6 DISCUSSION 

We have already verified that our main QPO frequencies agree with 
frequenci es obtained with an ind ependent, fiilly nonlinear numeri- 
cal code jCerda-Duran et ai]|20b7l) . for the same initial model. We 
caution, however, that we have not yet considered the crust-core 
interaction, different magnetic field topologies or the coupling to 
the exterior magneto sphere. These effec ts have to be taken into 
account and already ISotani et al.l ( 12007 j) , find that the observed 
QPOs could lead to constraints on the magnetic field topology. 
To complete the picture, a three-dimensional numerical simulation, 
that includes a proper coupling of the crust to the MHD interior 
and to the exterior magnetosphere will be required and our current 
results provide a good startin g point. Extensive de tails of our com- 
putations will be presented in lSotani et alj ( 12007 j) . 
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